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Abstract 



We study particle decay in de Sitter space-time as given by first order perturbation theory 
in a Lagrangian interacting quantum field theory. We study in detail the adiabatic limit of the 
perturbative amplitude and compute the "phase space" coefficient exactly in the case of two equal 
particles produced in the disintegration. We show that for fields with masses above a critical mass rric 
there is no such thing as particle stability, so that decays forbidden in flat space-time do occur here. 
The lifetime of such a particle also turns out to be independent of its velocity when that lifetime is 
comparable with de Sitter radius. Particles with mass lower than critical have a completely different 
behavior: the masses of their decay products must obey quantification rules, and their lifetime is 
CS| ' zero. 
> 

m 

iTi '. 1 Introduction 
rn 

04 ' Some important progress in the astronomical observations of the last ten years [1] [2] have led in a pro- 

gressively convincing way to the surprising conclusion that the recent universe is dominated by an almost 
spatially homogeneous exotic form of energy density to which there corresponds an effective negative 
pressure. Such negative pressure acts repulsively at large scales, opposing itself to the gravitational 
attraction. It has become customary to characterize such energy density by the term " dark" . 



. , The simplest and best known candidate for the "dark energy" is the cosmological constant. As of today, 

' the ACDM (Cold Dark Matter) model, which is obtained by adding a cosmological constant to the 

standard model, is the one which is in better agreement with the cosmological observations, the latter 
being progressively more precise. Recent data show that dark energy behaves as a cosmological constant 
within a few percent error. In addition, if the description provided by the ACDM model is correct, 
Fricdmann's equation shows that the remaining energy components must in the future progressively 
thin out and eventually vanish thus letting the cosmological constant term alone survive. 

In the above scenario the de Sitter geometry [31 H] , which is the homogeneous and isotropic solution 
of the vacuum Einstein equations with cosmological term, appears to take the double role of refer- 
ence geometry of the universe, namely the geometry of spacetime deprived of its matter and radiation 
content and of geometry that the universe approaches asymptotically. On the other hand, it seems 
reasonable to imagine that the presence of a small cosmological constant, while having a huge impact 
on our understanding of the universe as a whole, would not influence microphysics in its quantum as- 
pects. However this conclusion may have to be reassessed, because in the presence of a cosmological 
constant, however small, it is the notion of elementary particle itself which has to be reconsidered: 
indeed, the usual asymptotic theory is based on concepts which refer closely to the global structure 
of Minkowski spacetime and to its Fourier representation, and do not apply to the de Sitter universe 
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which is not asymptotically flat. Secondly, even if one may think that interactions between elementary 
particles happen in a "laboratory" so that "infinity" is a distance of the order of meters, our present 
understanding of perturbative quantum field theory is also based on global concepts; in particular, the 
calculation of perturbative amplitudes involves integrations over the whole spacetime manifold and it 
should be expected that different topological global structures result in different physical properties in 
the "small". 

The literature about de Sitter quantum field theory is very extensive, but there is no comparison with 
the understanding one has of Minkowskian field theory as regards both general and structural results 
as a well as its operative and computational possibilities. This second point is particularly doleful: 
calculations of perturbative amplitudes which in the Minkowskian case would be simple or even trivial 
become rapidly prohibitive or impossible in the case of de Sitter or anti de Sitter universes: this in spite 
of the fact that one is dealing with maximally symmetric manifolds which have invariance groups of the 
same dimension as the Poincare group. The technical, but also the physical, difference lies precisely in 
the above mentioned fact that much of the usual quantum field theory is based on concepts which are 
characteristic of the global structure of Minkowski spacetime and which do not persist in the presence 
of curvature, already in the presence of a mere cosmological constant, where Minkowskian spacetime is 
replaced by the de Sitter or by the anti de Sitter one. 

In this paper we give a full description of how to solve the problem of calculating the mean lifetime of 
unstable scalar particles on de Sitter spacetime at first order in perturbation theory. This interesting 
physical problem provides also an example of a concrete perturbative calculation in presence of the 
cosmological constant. The task already presents considerable mathematical difficulties. 

To our knowledge this calculation was first taken up by O. Nachtniann [5] in 1968. He showed, in a very 
special case, that while a Minkowskian particle can never decay into heavier products, a dS-particle can, 
although this effect is exponentially small in the dS-radius. 

The subject has acquired a greater physical interest with the advent of inflationary cosmology. In partic- 
ular, the idea that particle decays during the (quasi-)de Sitter phase may have important consequences 
on the physics of the early universe has been suggested recently [6l[71[8]. The mathematical and physical 
difficulties related to the lack of time-translation symmetry of the de Sitter universe, and more gener- 
ally of non-static cosmological backgrounds, have been tackled [6l[7l[8] by using the Schwinger-Keldysh 
formalism, which is suitable for studying certain aspects of the quantum dynamics of systems out of 
equilibrium. An important ingredient of this approach is the so called Dynamic Renormalization Group 
[9] which allows a kind of resummation of an infinite series of infrared diverging quantities. That method 
is however based on the introduction of a practical notion of lifetime of an unstable particle which is 
quite different from the definition commonly used in quantum physics. Also, the hard technical diffi- 
culties of the concrete calculation involved in solving a complicated integro-differential equation have 
only been faced in the favorably special conformal and minimally coupled massless cases although in 
principle the method can be used to deal with particles of generic mass [6j [71 18] . 

In this paper we perform a computation which is similar to the one outlined by Nachtmann and fol- 
lows the conventional quantum field theoretical perturbative approach for computing probability ampli- 
tudes. Our work gives significantly wider results w.r.t. ,5j, e.g. regarding the so-called adiabatic limit, 
complementary-series-particles, and explicit expressions of the relevant Kallen-Lehmann weights. On 
the other side comparing our result with those of [6, j8J is not easy because of the non standard (but 
interesting) definition of lifetime chosen in [H [7l [8] . 

These findings have been summarized in a recent short communication TO] . The results exhibit signif- 
icant differences compared to the Minkowski case, and decay processes which arc normally forbidden 
become possible and, vice-versa, processes that are normally possible are now forbidden. The maximal 
symmetry of the de Sitter universe implies the existence of a global square-mass operator, one of the 
two Casimir operators of the de Sitter group SOo{l,d) (see e.g. [11]); this quantity is conserved for de 
Sitter invariant field theories. However, in contrast with the Poincare group case, the tensor product 



Particle decay in the de Sitter Universe 



3 



of two unitary irreducible representations of masses mi and m2 decomposes into a direct integral of 
representations whose masses m do not satisfy the 'subadditivity condition' m > mi + ni2'- all repre- 
sentations of mass larger than a certain critical value (principal series) appear in the decomposition. 
This fact was shown in [5] for the two-dimensional case and will be established here in general. This 
means that the de Sitter symmetry does not prevent a particle with mass in the principal series from 
decaying into e.g. pairs of heavier particles. This phenomenon also implies that there can be nothing 
like a mass gap in that range. This is a major obstruction to attempts at constructing a de Sitter 
S-matrix; the Minkowskian asymptotic theory makes essential use of an isolated point in the spectrum 
of the mass operator, and this will generally not occur in the de Sitter case. We will also show that the 
tensor product of two representations of sufficiently small mass below the critical value (complementary 
series) contains an additional finite sum of discrete terms in the complementary series itself (at most 
one term in dimension 4). This implies a form of particle stability, but the new phenomenon is that a 
particle of this kind cannot disintegrate unless the masses of the decay products have certain quantized 
values. Stability for the same range of masses has also been recently found [8] in a completely different 
context. Other remarks about the physical meaning and applicability of our results will be presented in 
the concluding section. 



1.1 Notation 

We denote C+ = — C_ the open upper complex half-plane. Let A = C \ [— 1, 1], Ai C\{~oo, 1]. The 
function log is defined as holomorphic on C\ (— oo, 0] and real on (0, +oo) and C,^-^ C,^ &s exp(/ilog(C)). 
It is entire in ^x. If C e Ai and p > 0, then {pC,Y = p^C^- If C e C+ and s e C_, then (sC)'' = s^'C^. 
We define z ^ (z^ — 1)^/^ as holomorphic on A and asymptotic to z at large \z\. It is Herglotz, negative 
on (— oo, —1) and positive on (1, -|-oo). 



2 Free fields in Minkowski and de Sitter spacetimes 

In this section we give a short summary of the theory of free and generalized free quantum fields on 
de Sitter spacetime. Since there are infinitely many inequivalent representations of the field algebra, a 
(mathematical) choice has to be made on physical grounds. Ours is based on the analyticity properties 
of the vacuum expectation values: see the condition (W2) below. In the Minkowski space, this is 
equivalent to the positivity of the energy. In the de Sitter case, it admits a thermal interpretation 
(T2, 13, 14l[l5]). The reader can find in [T3l[T4l[T5] a general approach to to de Sitter QFT based on 
such analytic properties. It includes the so called Bunch-Davies, also called Euclidean vacuum of de 
Sitter scalar Klein-Gordon fields as a basic example. 

The real (resp. complex) d-dimensional Minkowski spacetime Md (resp. Mjf'') is R'' (resp. C^) equipped 
with the Lorentzian inner product 

/ _ /O 1/1 d—l fd—l Q fO ^ ^/ {o T\ 

3C ' 3j — ijC 00 00 ■ ■ < 00 00 — 00 OC 00 ' 00 ( ^ < .L ) 

w.r.t. an arbitrarily chosen Lorentz frame {e^,a = 0, . . . , d — 1} . When no ambiguity arises, =^ 
X ■ X. The real (resp. complex) de Sitter spacetime Xd (resp. xjf ^) with radius R > are the 
hyperboloids 

Xd^{x€Md+i : x-x + R'^ = 0}, X^/'> ^ {x £ M^/_l^ : x-x + R'^ = 0}, (2.2) 

equipped with the pseudo-riemannian metric induced by ()2.ip . L\_{d) = 5*0(1, d—l; R) is the 
connected Lorentz group acting on Md, and i+(C; d) — 5*0(1, d—l; C) is the connected complex 
Lorentz group acting on M^^-*. The connected group of displacements on Xd (resp. xj^'^^) is L]^_{d + 1) 
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(resp. i+(C; d+ 1)), sometimes denoted Go (resp. Gq^^). These groups act transitively. Note that 
our definition of etc. arbitrarily selects a particular orthonormal basis (eo, ea-i) in Md or 

(eo, . . . , Cd) in Md+i- These particular Lorentz frames will be useful in the sequel. In the future 
and past open cones V± and the future and past light-cones C± are given by 

V+ = {x e Md : X ■ X > 0, x° > 0} = -V- , 

C+ = {x G Md : X ■ X = 0, x°>0} = -G_ . (2.3) 
The future and past tubes in the complex Minkowski spacetime M^"^ are given by: 

T± = + iV± . (2.4) 

(c) (c) 

The future and past tuboids in are the intersections of the future and past tubes in M^_^-^ with the 
complex de Sitter manifold xj^'^: 

T± = T± n X^"^ . (2.5) 

We will use the letter X to denote cither Md or Xd when the same discussion applies to both, A"*^^^ 
denoting the complexified object, dx will denote the standard invariant measure on X , i.e. using the 
frame (eg, . . . , e„), dx = dx^ . . . dx"^^^ in the case of Md, and dx — 2(5(.T^ + i?^) dx" . . . dx"^ for Xd- 

A (neutral scalar) generalized free field on A" is entirely specified by its 2-point function. This is a 
tempered distribution W on X x X (we denote W{x, x') = W{x', x)), which we require to have the 
following properties: 

(Wl) Hermiticity: 



W{x, x') = yV{x', x) . (2.6) 

(W2) Analyticity and invariance: there is a function w of one complex variable, holomorphic in the cut 
plane C \ R+, with tempered behavior at infinity and at the boundaries, such that, in the sense 
of tempered distributions, 

yV{x, x') = lim w{{z - z'f), (2.7) 

z e T_, z' e T+ 

hence 

W'{x, x') = W{x',x) = lim w{{z-z'f). (2.8) 

z e T+, z' G T_ 

For complex z, z' G X^""^ such that {z - z'f e C \ R+ we will denote W{z, z') = w{{z - z'f). 
Note that this implies 

W{z, z') = W{z', z) = W{-z', -z), (2.9) 

and 

W{x, x') = W{-x', -x) . (2.10) 

(Wl) and (W2) also imply 

W{z, z') = W{z, z'). (2.11) 

(W3) Positivity: For every / e S{X), 



{yV,f<»f)= f{x)W{x, x')f{x')dxdx' >0 . (2.12) 

JxxX 
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Conversely, given W and w having these properties, after having identified the kernel A/i — {f E S{X) : 
(W, f ® J) = 0} we can construct a Hilbert space Ti by completing S{X)/Mi equipped with the scalar 
product (/, g) = (W, f ® g), and then exponentiate J-'i into a Fock space 

oo 

= jr^^ jr^ ^ C, = 5J^f " for n > 1 . (2.13) 

n=0 

The vacuum is the unit vector 1 G — C. There is a continuous unitary representation U of the 
Poincare or de Sitter group acting on and preserving the J^m with Ufl = fl. The generalized free field 
(j) is defined on a dense domain in !F and {il, (j){x)(j){x')il) — W{x, x'). 

As a result of the analyticity property (W2) , the Wick powers of a generalized free field are well-defined 
local fields operating in the same Fock space. Their vacuum expectation values are obtained by the 
standard Wick formulae as sums of products of W. 

We note that a function W on Xd x Xd possessing the properties (Wl) and (W2) automatically extends 
(through (|2.7p ) to a function with the same properties on Md+i x Md+i, so that a generalized free field 
on Xd has an extension as a generalized free field on Md+i- However the extension of W need not satisfy 
(W3) on Md+i X Md+i even if it does on Xd x Xd- 

A free field (j) of mass m > on A:" is a generalized free field such that W is a solution of the Klein-Gordon 
equation with mass m in both arguments, and is normalized so as to obey the canonical commutation 
relations. In that case W is uniquely determined by m and will be denoted Wm- In the Minkowskian 
case, the representation U\J-i is irreducible and equivalent to the representation [m, 0] of the Poincare 
group. 

As usual, the representation U provides a representation of the Lie algebra of the (Poincare or de Sitter) 
group and its envelopping algebra by self-adjoint (or ix self-adjoint) operators on T. In particular the 
square-mass operator is given by = P^P^ in the Minkowskian case, and by = M'^" M^,y/2R'^ 
in the de Sitter case. In both cases, M^^* = m^Vl/ for every £ T\. See e.g. [Tl] . 

2.1 Special features of free fields in de Sitter space-time 

In the de Sitter case the mass m can be related to a dimensionless parameter v as follows 

n?R^^('^-^S + y\ (2.14) 



2 

-| 1/2 

= ±i?(m2-m^)i/2, m, = ^. (2.15) 



In this case, if no ambiguity arises, we shall often denote = >V_^ to mean Wm, and similarly 



and Wl,. Explicitly, if z, z' €E xj^'^\ [z — z'Y ^ R-+, and hence Q ~ z ■ z' jP? does not belong to the real 



interval — oo, — 1]), 



wA-^.z') . (c^-i)-^p-fy(c) (2.16) 



(4^)l^d-2r(|) V 2 '2 ' 2 ' 2 

Since T{c)~^ F{a, b ; c ; z) is entire in a, b, and c the rhs of (|2.17p is meromorphic in ly with simple poles 
at ^ = zti{{d— l)/2 -I- n) , n > an integer. In other words Wi, {z, z') extends to a holomorphic fmiction 
of z and z' in the domain {j/ e C, z G z' e X^'^^ : i/ ^ ±i((d - l)/2 + Z+), (z-z')2^R+}. 

However w^, possesses the positivity property (W3) (see (I2.12p ) only if either 
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(1) V is real, i.e. m > rric — {d ^ l)/2R. In this case U\J-i is an irreducible unitary representation of 

the "principal series" . 

or 

(2) I' is pure imaginary with iiy E (— (rf— l)/2, (d — l)/2), i.e. < m < (d— 1)/2R. In this case U\J-i 

is an irreducible unitary representation of the "complementary series" . 

We shall need a small part of the harmonic analysis on the de Sitter space-time as developed in |14| . 
If z e T± C X^f^ and ^ S C+ \ {0} C M^+i, then ±Im(z ■ > 0, so that (z • is well-defined and 
holomorphic in (z, A) in (7+U71) x C. The role of plane waves on Xd is played by the distributions 



i^t{x,0^ Hm ({x±ty)-0^^i^f{x,0- (2.18) 

An important formula expressing the de Sitter case two-point W^, as a Fourier superposition of plane- 
waves is the following (see [14]): 

W.{z, z') = Rca,. f {z ■ O-^^'-" • z')'"^"" «(C), (2.19) 

where zi e 71 , Z2 G 7+ , and 

C'^- = ^^irri^ ■ (2-20) 

In ((27T9)) . 7 is a (d - l)-cycle in C+ \ {0} homologous to the sphere Sq = C+ n : = 1}- The 
(d — l)-form a is given, in the standard coordinates, by 

d 

a = (e°)-i X!(~l)'^'^'' • ■ ■ '^^^ • • ■ ^^'^ ■ (2-21) 

If a smooth function / on C+ \ {0} is homogeneous of degree (1 — d), the form fa is closed, so that the 
linear functional 

/ ^ W) - / /(0«(e) (2.22) 

is independent of 7. This implies that it is Lorentz- invariant. We often denote dfx-y the measure defined 
on 7 by the restriction of a. In particular the restriction of a to the {d — l)-sphere Sq is the standard 
volume form on that sphere, normalized by J^^ dfiSgiO — 2n'^^^ /r{d/2). It is possible to take the limit 
of (|2.19p . in the sense of distributions, when zi and Z2 tend to the reals: 

W,(x, x')^Rcd,. f V^:.-!^, (a;,OV^+.-i_, (x',e)rf/z^(0 . (2.23) 

Jj 2 '^'■^ 2 

Comparing ([2171) with ([219]) and ([2:201) gives 

,.-^+..,, w e-^2^ d-1 . d 1-C 

(^•0 ^ (e-2) '^^'^(^^ " ;^^^TF(|y ' 2 ' ^ 

(2.24) 

Both sides of this equation are holomorphic in zi, Z2, in the domain x T^i- x C, hence the equation 
(|2.24p holds in this domain. 

Remark 2.1 If T is a homogeneous distribution of degree /3 on C+ \ {0}, it can be restricted to any C°° 
submanifold of dimension d — 1 which is transversal to the generators of C+, in particular to hyperplanar 
sections such as Sq = {£, e C+ : = 1} and Vq = {£, £ C+ : + ^'^ = 1}. If 7 is of this type and 
compact, T{£^)a{£,) is well-defined and, if /3 = 1 — d, it is independent of 7. 
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Remark 2.2 For any complex a, (z, ^) i-^ (z • is C°° in ^ and holomorphic in z on 7± x (C+ \ {0}) 
and it is an entire function of a. For each ^ it has a hmit in the sense of tempered distributions on Xd 
as z tends to the reals, and this has been denoted il)^{x, ^). It is an entire function in a. Furthermore 
its invariance under Go implies that, if g S{C+ \ {0}), /^^ "^a (^j OviC) '^^ is i'^ ^- Indeed any 
small displacement of x can be effected by a group transformation close to the identity, which can be 
transferred to ^ and thence to (p. In the same way, C) is C°° in ^ (as well as homogeneous) 

when integrated with a smooth test-function in x. This explains the meaning of formulae such as 
(I2.23p . Note that the integral in this formula is entire in v. For similar reasons, for any e iS(Xrf), 
ipix) yV,y(x, x') dx is C°° in x' and meromorphic in v. 

2.2 More features common to Minkowski and de Sitter space-time 

An important formula, which holds in Minkowski as well as in de Sitter space-time (but in this case 
only if m, m' > mc), is the projector identity: 

[ Wmiz, x)Wm'{x, y)dx^Ci{m, d)S{m^ - m'^)Wm{z, y). (2.25) 
Jx 

Here 

Ci(m, d) — 2tt for Minkowski space-time, (2.26) 
Ci(m, d) = Co (i') = 27r| coth(7riy) I for de Sitter space-time. (2.27) 

The proof of the above identity is trivial in the Minkowskian case. For the de Sitter case it will be 
provided in Appendix [Pl Note that Co{mR) tends to 27r as i? ^ +oo for a fixed m > 0. 

The K alien- Lehmann decomposition theorem exists in both and Xd- In the case of Md, (see jl6) . 
p. 360), it asserts that, for every W having the properties (Wl) and (W2) there is a tempered p such 
that 

W{z, z')= I p{m^)W^{z, z')dm^ . (2.28) 

If W satisfies (W3), then p is a tempered positive measure. The same holds in the dS case provided 
W satisfies some decrease property. In this case, the integral runs on masses of the principal series, i.e. 
m > TOc = (d — 1)/2R. For proofs and details, see [HI [17]. In particular if mj > and, in the dS case, 
ruj > mc for 1 < J < iV, 

[|W™^.(x, x')^ / p{a^; mi,..., mjv)Wa(x, x')da^. (2.29) 

j=l J a>b 

Here b = mc in the de Sitter case, b = J^j ™j in the Minkowski case. 

3 Particle decays: general formalism 

There is at the moment nothing like the Haag-Ruelle asymptotic theory (HRT) (see [lH [TH [16]) for 
the de Sitter universe. Indeed all the ingredients of that theory are missing in the de Sitter case. For 
example, as it will be shown in this paper, even in a free field theory of mass m > mc, the mass m is not 
an isolated point in the mass spectrum. Moreover the solutions of the Klein-Gordon equation do not have 
the kind of localization at infinity which plays an essential role in the HRT. The concept of a particle is 
therefore not obvious in de Sitter space-time, except for localized observations. Here we adopt Wigner's 
point of view: a one-particle vector state is a state belonging to an invariant subspace of the Hilbert space 
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in which the representation of the invariance group reduces to an irreducible representation. In the dS 
case, we also require that this irreducible representation belong to the principal or complementary series, 
i.e. it should be equivalent to one of the representations which occur in the T\ of a free field. 

We shall study the decay of a particle using first-order perturbation theory. The initial framework 
and calculations are the same for the Minkowski and de Sitter cases: its ingredients are the projector 
identity and the Kallen-Lehmann representation. (It can also be extended to the Minkowskian thermal 
case ([20]) although there is no Kallen-Lehmann representation there). Let 

00, 01, ■ ■ • , 0Ar (3.1) 

be 1 -|- A'^ independent free scalar fields with masses toq > 0, m\ > 0, . . . , m^v > 0, acting in a common 
Fock space 7i, the tensor product of the individual Fock spaces for the (/)fc : 



w-(g)^''•■^ (3.2) 

fe=0 

(SI. 4>,{x) My)^) = S,k W™, (x, y) . (3.3) 

We denote 

w.o,....«=-^r^---^-^if • (3-4) 

This is the subspace of states in Ti. containing jk fc-particles. ^/joi-.-jN denotes the hermitian projector 
onto this subspace. We now switch on an interaction term 

-fg{x) C{x) dx, C{x) = : Mx)<Pi{xf^ . . . 0w(x)'" : . (3.5) 

X 

Here the qj are non-negative integers, and we denote q\ — JljLi'Zi'- 7 is a small constant, gr is a 
smooth, rapidly decreasing function over X . In the end, g should be made to tend to 1 (adiabatic limit). 
According to perturbation theory, the transition amplitude between two normalized states '00 and ipi 
in TL is given by (?/;o, S{'yg)ipi), where S{'^g) is the formal series in 75 

'S'(75) = — / g{xi)dxi---g{xn)dxnT{C{xi)...C{xn)) (3.6) 

In (|3.6p . T{£{xi) . . .£(a;„)) denotes the (renormafizcd) time-ordered product of £(xi), . . .£(a;„). In the 
first order in jg, the transition amplitude between two orthogonal states ipo and ipi is 



ii^o, iTi{'-fg)ipi), ^1(73)=/ -/g{x)C{x)dx . (3.7) 

We take 



X 



Ipo = j k{x)(t)o{x)Q.dx , (3.8) 

ipi = / /i(a;ii, . . . . . . , SATi, . . . ,a;Arq„) : ]^ J]^ 0j(xjfe)(ia;j7c : , (3.9) 

j=i fc=i 

where /o and /i are smooth rapidly decreasing functions. The states of the form (|3.8p generate 7ii,o...,o 
and the states of the form (|3.9p generate Ho,gi,...,gjv ■ The probability of transition from i/jq to any state 
in Ho,gi,...,g„ is: 

^ _ (tAo, ri(7g)i;o,,„...,g,Ti(7g)* V>o) ^ qlj^ 

(00, V'o) (-00, -00) 



foi^) fa{y) 9{u) g{v) X 
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N 



X Wm^ix, u) |n>V™^.(w, vy^jWmoiv, y) dx du dv dy . (3.10) 

From now on, we suppose, in the dS case, that rrik > rric, < k < N , i.e all particles belong to the 
principal series. We may then replace the central two-point function in u and v by its Kallen-Lehmann 
decomposition: 

N 

JJWm^-(u, v^' = / p(a^;mi,...,TOi,...,mAr,...,mAr)Wa(u, v)da'^ . (3.11) 
Here rrij occurs Qj times as an argument of p. This gives 



(*o,*o) 



foix) foiy) giu) g{v) p{a ; mi, . . . , wi, . . . , mw, . . . , mAr) x 

X VVmo(x, u)Wa{u, v)Wr,io{v, y) dx du dv dy do^ . (3.12) 

The next step would be the so-called adiabatic limit, and should consist in letting the cut-off g tend to 
1 in this formula. It is however easier to set first only one of the g's equal to 1, say g{u) = 1 in (|3.12p . 
It then becomes possible to perform the integration over u by using the projector identity (|2.25p and 
we find for the transition probability: 

r = Li{fo, g) X q\p{ml]mi,. . . ,mi, . . . .ttzat, . . . .mjsi) , (3.13) 

where 



^^^j^ ^ I'^Ciijno, d) J g{v)fo{x)Wmo{x, v)Wmo{v, y)fo{y)dx dy dv 

J foix)Wmo{x, y)fo{y)dxdy 
This formula exhibits an interesting factorization: the first factor depends only on the wavepacket /o, 
the mass mo of the incoming particle and the switching-off factor 7^(7; the adiabatic limit still remains 
to be done there; the second factor contains all the information about the decay products. 

If we now attempt to set g{v) = 1 in (|3.14p and to integrate over v using again (|2.25p . the result is 
proportional to 6{mQ — mg), i.e. the integral diverges. This difficulty was resolved in the 1930's by 
aiming at the average transition probability per unit time (see e.g. j21], pp. 60-62). We first review the 
well-known Minkowski case, in a form which can serve as a model for the de Sitter case. In fact even this 
famous old case deserves some re-examination on its own right and it is possible, in this case, to allow 
the two g in (|3.12[) to tend to 1 simultaneously, or even at different rates. This is done in Appendix 
lAl It is found that, if both g are taken as in (|4.3p . the result is the same as found above. But this is 
not necessarily the case for other g. Nevertheless the procedure announced above (i.e. setting the first 
g in (|3.12p be equal to 1, then discussing the time average of the limit as the second g tends to 1) will 
be used in the de Sitter case, since it gives good results in the Minkowski case, and since calculations 
in the dS case would become much more difficult otherwise. Note that in the de Sitter case (|3.13p and 
p.l4p are applicable only when mp > rric and the range of integration over in (|3.12p contains only 
values > ml {rric = {d— l)/2/R). In the case of the decay into two particles of mass mi, it will be 
seen below that this includes the case mi > mc, but also the case nic > mi > mc\/3/2. 



4 Minkowski case 

4.1 Adiabatic limit: the Fermi golden rule 

The simplicity of the Minkowskian case arises from being able to use of the Fourier representations: 

fo{x) = J e-^P^fo{p)dp, g{x) = J e-'P^g{p)dp, w^{x, y) = {2^f-^ j e'P^y-^U{p' - my{p')dp. 

(4.1) 
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Then the factor in p.l4p becomes 



^ (27r)^7^ / /o(p) Sjp^ - mg)g(pP)/o(g) Sjg^ - ml)e{q^)g{p - q) dpdq 

We now speciahze the cut-off g to depend only on the time coordinate of the chosen frame g{v) = 
hiv'^) = h(t), i.e. we think of the interaction as smoothly switched on and then turned off. The Fourier 
representation is then g{p) — h{p^) d{p) and cq. (|4.2p becomes 

_ i2n)wm m'r' imi' s{p' ^ m^)g(/) dp 

J\foipWSip^-ml)e{pO)dp 

If we choose for g the indicator function of a time-slice of thickness T, i.e. h{t) = 9{t + T/2) 6{T/2 — t), 
h{0) = T/2n, we get 

. _ y / (27r)7^ /(2/)-^ \Mp)? Sjp' - mDejp^) dp \ 
°' V J\fo{pWS{p^^ml)9{pO)dp J ' 

Therefore, as noted above, removing the cut-off produces infinity. However, according to the Fermi 
golden rule, what is physically meaningful is not the amplitude but the amplitude per unit time. There- 
fore, dividing this by T and taking the limit as T — s- oo (a particularly trivial operation in this case) we 
finally get the following expression for the transition probability per unit time: 

1 _{2n)j^ J{2pO)-^\Up)\^ 6{p^-ml)e{pO)dp , ^ ^ 



^(/o) J\fo{pWS{p^-ml)e{p")dp 



ql p(mo; mi, . . . , wi, . . . , ttin, ■ ■ ■ , mN)- (4.5) 



The reciprocal of this expression is the lifetime of the 0-particle in the state /q. The dependence on 
the wavepacket /o is a crucial feature of the special relativistic Minkowski case as it will be readily 
recognized. For instance to compute the lifetime tq of a particle at rest in the chosen frame we may let 
\k{p)? tend to 5{p), e.g. by taking 

7o(p)=e(i-'^)/V(p/e), h{x)^2^5{x^)e^''-''^'\{ex), e > 0, (4.6) 
with ip e SiBJ^-'^), and letting e ^ 0. Then tends to 

— = — ^ q\ p{ml; mi, . . . , mi, . . . , niN, ■ ■ ■ , mN) ■ (4.7) 

To mo 

We may act with a Lorentz boost on the same particle by replacing in (|4.5p the wavepacket /o by 
fgix) = fo{A^^x), A e the amplitude is modified as follows 

1 _ (2^)72 /(2(Ap)0)-i |/o(p)|2 S{p^ - ml)e{p°) dp ^ . , 



^ifo) J\fo{pWS{p^-ml)9{pO)dp 



(7!p(mo;mi,...,mi,...,mAr,...,mAr) . (4.8) 



If again |/o(p)P ^{p)i tJ^c final result is the expression in (|4.7p multiplied by l/(A)oo. If A = 
exp(sAfio), i.e. the particle is moving with velocity v — ths, (A)oo = coshs = (1 — v^)^^/^ gives the 
usual correction to the lifetime: 



ro/Vl-v2. (4.9) 



Remark 4.1 It is worthwhile to stress once more that this effect, which expresses the behavior of the 
life-time of a moving particle in special relativity, crucially depends on the peculiar way in which the 
wavepacket enters in the transition amplitude per unit time (|4.5p . 



Particle decay in the de Sitter Universe 



11 



4.2 Kallen-Lehmann weights 

The weight p can be expHcitly computed only in the case of one particle decaying into two particles. 
For a particle of mass toq > decaying into two identical particles of mass mi > 0, i.e. the case N = 1, 
gi = 2, the well-known formula is 

d-3 

p{ml; rm,mi) = '-j— ^(m^ - 4m?), (4.10) 

{in)^2d-2r(^) mo 



and (14.71) becomes 



d-3 

— = (lifetime of mo)" = b^-^ 0{ml - 4m?). (4.11) 

^" (4^)^2d-3r(^) m2 



For d = 4 this is 

,2 ,1 ™ 2 



To OTT mg 

in agreement with the computation in e.g. |21j . 



1 7^ (m? — 4m?) ^ , ., , , 

- - ^ ° 2 g(m^-4mf), (4.12) 
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5 de Sitter case 

5.1 Adiabatic limit in the de Sitter case 

The discussion of the adiabatic hmit is more comphcated in the de Sitter case. Taking the adiabatic 
hmit is of course technically much more involved than in the Minkowski case (and we will relegate all the 
technical details to the appendices). But the really intricate and maybe perplexing issue is the physical 
interpretation of the whole procedure and, even more, of the somewhat surprising results. 

Having in mind the Minkowskian case that we have just discussed, the first question that should be 
asked is what is "time" in the de Sitter universe and what does it means that an interaction lasts for a 
certain time. In the Minkowski case we have the solid foundation of special relativity and a privileged 
class of frames, the inertial frames, each of them having an inherent precise notion of time. 

In the de Sitter case (and the situation is even worse in a general curved spacetime) we have no such 
thing. Instead we have many possible coordinate systems, that may or may not cover the whole manifold, 
and many possible choices of temporal coordinates that have no special relation to each other. 

For example, the de Sitter universe is the only known spacetime manifold admitting three different 
inequivalent choices of cosmic time so that the de Sitter metric takes the appearance of a, respectively, 
closed, flat, or open Friedmann-Robertson- Walker universe. But there are also other possibilities. The 
choice of time coordinate made in 1917 by de Sitter in his original papers [3l |4] describes a wedge-like 
region of the de Sitter manifold as a static spacetime with bifurcate Killing horizons [22] . 

We choose to proceed heuristically in analogy with the Minkowskian case. Concretely, we will work out 
the adiabatic limit using two of the three possible cosmological coordinate systems, namely the closed 
and the flat systems. Starting again from eq. p.l4p we take the cutoff g appearing in there as the 
indicator (or characteristic) function of some "cosmic time-slice" of thickness T w.r.t. to the relevant 
choice of cosmic time. 

We will see that in both the closed and flat case the amplitude diverges linearly in T precisely as in the 
Minkowskian case. Therefore, to extract a finite limit we are entitled (and have no other choice than) 
to use the Fermi golden rule and compute in the above two frames the probability per unit time by 
dividing by T\ there is at this point a small difference w.r.t. the flat case: the amplitude per unit time 
at finite T depends on T. However, letting T — > oo gives a well-defined limit which exhibits a much 
more disturbing difference with the Minkowskian case. 



Closed FRW model: The relevant coordinate system is the following: 

, , J = R sinh(i/i?), , , 

^^ '"^ \ ^ Rcosh{t/R) i-u), u£S'^-\ ^ ' 

(the minus sign at rhs is for further convenience) . In this coordinate system the constant time slices are 
hyperspheres. These coordinates have the advantage to globally cover the de Sitter manifold; they gives 
to the metric the form of a closed FRW model with scale factor a{t) ~ cosh{t/R): 

ds^ = dt^ - R^cosh^{t/R)da'^{u) (5.2) 
= dt^ ~R^cosh^{t/R){dx^ +siii^x{dO^ +sm^Od(l)'^)) (rf = 4). (5.3) 

In dOl) da^{u) is the square line element on S'^ ^ at u, and (|5.3p includes its expression in Euler angles. 
In these coordinates, we choose 



g{x) = grix) = 9{t + T/2) 0(T/2 - t) 



(5.4) 
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Figure 1: Time-slices of the de Sitter spacetime in the closed and in the flat coordinate systems. 



Flat FRW model: These are the coordinates currently used in the context of inflationary models. 
Hypersurfaces of constant time are flat: 



R sinh + 2l? ^ ^ 
e « y-' , ( 1 < i < c? - 



R cosh ^ - 



1) , y e R'' 



(5.5) 



idyl 



In these coordinates we choose 



-lix) = grix) = 9{t + T/2) 0{T/2 - t) 



(5.6) 



(5.7) 



But the coordinates (|5.5p only cover one half of Xd, the region where + x'' > 0, and the adiabatic 
limit will have to include the contribution of the other half not covered by the coordinate system. 

It turns out that the limit 



1 ( y Li{fo, .g)\ , 2 

- = 1™ 7^ X 'l-P{mo;mi, ...,mi 

T \T^oc 1 



,mN, 



,mAr) 



(5.8) 



exists and is the same for both kinds of slices. The calculations are tedious and not quite straightforward, 
and will be given in Appendices [B] and For the spherical slices of the closed FRW system, only the 
calculations for d = 2, 3, 4 have been carried out. The method for the flat FRW coordinates works for 
all d. 



The inverse lifetime that results is 

^l,qi,...,qN 



(Jo) = n ^ g!p(™o;™i' ■ • ■ . . . ,mjv. 



,mN) 



(5.9) 



where we denoted n = R{mQ — m^)^/^. Note the similarity of this formula with (|4.7p and in fact the first 
factor in (|5.9p tends to the corresponding factor in (|4.7p when i? ^ oo at fixed mo . However there is a 
most striking difference with (|4.5p : the rhs of (|5.9p does not depend on /q, the initial wave function of 
the decaying particle. Therefore in particular the lifetime of a particle does not depend on its velocity. 
We will comment on this feature, at first sight embarrassing, in the conclusions. 
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6 Kallen-Lehmann weights 

As in the flat case, an explicit computation of the KaUen-Lehmann weight is only possible for decays 
of one particle into two. Here the discussion will be restricted to the case of a particle of mass mo > 
decaying into two particles of equal masses nii = > and we suppose at the beginning mi > rric. 
The more difficult case mi ^ m2 will be treated in a paper in preparation [23j . We shall find an explicit 
p{a?]mi, mi) such that 

Wm{z, z'f = / p(a2;mi, mi)Wa{z, z')da^ . (6.1) 



We change to variables k= [a? — {d—lY / , v = [m\R^ — {d— 1)-^/4]^/^, and (by abuse of notation) 
seek a function p{k] v) (mostly abbreviated as p{n)) such that 



VF^(z, z'f= I 2kp{k-v, iy)W^{z, z') dn . (6.2) 
By ()2.16p . this is equivalent to 

C'l, (x2 - ir'^PZIlJxf ^ I 2C',^^kp{k) PZIlJx) dn , (6.3) 



with 

T + ii^)T -iiy 



2(27r)ii?'^-2 

The generalized Mehler-Fock theorem ([H], p. 398) asserts that 



/■OO 

9{^) = / P\+,,{x)f{K)d^ 
JO ^ 



'0 

/(«) = ^sinh(^K)r Q -d + m') r Q -a-zK ) / PZi+,,{x)g{x)dx . (6.5) 



Therefore (|6.3p implies 



kp{k) = ^^l" ^sinh(7rAt) T ^— y— + T - ii^^ x hdin, v), (6.6) 



J I \ dot / ,2 1 N — — 

hd{K., V, v) ^ [X ~l) -1 



P ^Z ix) P (x)dx . (6.7) 



It is possible to obtain an explicit expression of /id(K, v) by using Mellin transform techniques (see 
|25j ) and a lemma of Barnes (see [H]). Recall that if e 25((0, oo)), its Mellin transform (p is given 

by 



oo 



^(s)= / C-V(C)rfC, (6.8) 
Jo 

It is entire in s = cr + it, decreasing faster than any negative power of r for fixed cr, and 

^(C) = — / r^^(s)rfs VaeR, (6.9) 

V'(C)=¥'(1/C) ^ ^(s)=^(-s). (6.10) 

If ifi, are in X>((0, ex:)), 

■^(0=/ <^l(C/")<P2(u) — <^(s) = <^i(s)(^2(s) . (6.11) 
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In particular (Mellin-Planchcrcl identity) 



(filiu) (p2{u) 



du 



1 

2iTT 



ifi(-s)if2{s)ds . 



(6.12) 



These properties can be extended to other functions and generahzed functions (see |25|), and, in many 
interesting cases, ahhough the MeUin transforms are no longer entire, the above formulae survive pro- 
vided the integration in (|6.9p or (|6.12p is performed on a suitable contour. 

we find 



By the change of variable x = y/1 + ( in 



hdi^, v,v)= Gi(C)G2(C) 



with 



and 



Gi(c) = ^'}^^\ G.(c) = c-^\pi;{^/^of 



a — + in, P — - 

The Mellin transforms of Gi and G2 are known (see [2S], 17(1) p. 257, and 28(1) p. 263.) 



Gi{s) 
provid 

G2{s) 



2^-1 



2' ' 2 ' 2 



r(i + ^^)r(i^^ 

provided Re/i < 2Res < min{2 + Rea, 1 — Rea}, 
1 



(6.13) 

(6.14) 
(6.15) 



7r2r(l + /3-/i)r(-/3-M) 



■ r 



1 - f + s, /3 - H - s, -/3 - 1 - ^ - s, -i±ii - s 



provided Re/i < Re(s + 1 — /i/2) < min{— Re/3, i}, 



(6.16) 

5 

(6.17) 



Remark 6.1 We have actually checked the above formulae using the methods described in [25]. How- 
ever some other formulae appearing in that extremely useful reference have misprints. 



By the Mcllin-Plancherel theorem 



hd{K, v, !/) = — 
ZlTT 



cr+zoo 



Gi(s)G2(-s)ds , 



and the preceding formulae give 

Gi{s)G2{-s) = 



2"^ I IK d-\-l iK d— 1 



iv. ^ - iiy] 



It is now possible to use Barnes' Second Lemma p. 112) 
Lemma 6.1 (Barnes) 



1 
2i7r 

= r 



ds 



provided 



fll + S, 02 + s, 03 + s, 61 — s, 62 — s 
c + s 

ai+61, 02 + 61, 03 + 61, ai + 62, 02 + 62, 03 + 62 
c — oi, c — 02, c — 03 

oi + 02 + 03 + 61 + 62 - c 



(6.18) 



(6.19) 



(6.20) 



(6.21) 



and that the contour 0/ integration in i6.20\) separates the increasing and decreasing series of poles. 
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As a function of s, Gi{s)G2{—s), as given by (|6.19p . is, up to a factor, of the form of the integrand of 
([OOl) if we take 



d-4 d-4 d-4 3 iK 3 iK 3d -6 

This choice satisfies the condition (|6.2ip . Therefore 

1 



(6.22) 



hd{K, V, v) 



rf— 1 \ in \ ■. d— 1 1 ■,. d— 1 I in d— 1 mi ■,. d~ 1 -,. d— 1 in 

p 4 ' 2 ' 4 ' 2 4 ^ 2 ' 4 2 ^ 4 2 ''^'4 2 

From here on, we will use the notation fi = {d — l)/4. 

We can recast the above expression for hd{K, v, v) by using Legendre's duplication formula: 



(6.23) 



1 



23-d/27r3/2r(2^ + iK)r(2/i - iK)T{2fi + iiy)^T{2fi - iiy)^r{2fi) 
in \ ^ f in \„/ ikX^ 

r ( + - + zz/j r u + y - *H r U + y ) x 



(6.24) 



In this form the formula is a special case of the formula for two unequal masses which will appear in 
[25] . We note that in the derivation of (|6.23p or (|6.24p with defined in (|6.7p , d is not restricted to be 
an integer. These formulae hold wherever both sides are defined. Eqs. (|6.23p and (|6.6p give 



kp{k; V, v) = 



K sinh(7rK) 



2d+^n^Rd-^r (1 + ^ + f ) r (i + /i - f ) r(2^) 



iK \ / iK \ / iK 



(6.25) 



or, using Ksinh(7rK) = 7r[T(iK)T{—iK)] ^, 
K p(k] V, v) = 



1 



2^+2,, 4- ^d-2 t{ik)T(^ik)T (i + ^ + f ) r (i + ^ - f ) r(2/i) 



e=±l 



' = ±1 



This obviously extends to an even analytic function of k, hence 

w^{z, z')'^ — / kp{k]v, v)wk{z, z')dK 



(6.26) 



(6.27) 



Moreover, for real v and k ^ 0, k p{k] v, v) is strictly positive. This shows that, in the presence of 
a suitable interaction term (see (j3.5p ). any "principal" particle can decay into any pair of equal-mass 
"principal" particles. 



6.1 Minkowskian limit 

Setting K — MR > and v = mR > in (|6.25p gives 

sinh(7rMi?) 



p{MR;mR, mR) 



2d+2^4^i?d-2r(d^) 



nLir(a;j +iRuj) 



T{xi + iRui) 



(6.28) 
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ith 



d-1 d+1 



M M MM , , 

= — +m, U2 = — -m, -U3 = — , U4 = — . (6.29) 

Recall Stirling's formula ([27], p. 47): 

r(z) = (27r)3e-^+(^-3)i°g^ (l + a^;j-i+a2^"^ + 0(2"^)), ai = 1/12, aa = 1/288 , (6.30) 

valid for z ^ R_. By a straightforward calculation it follows that if z ~ x + iy and \x\ remains bounded 
while |?/| — > +0O, 



IT.. , ■ m2 o -Trial I \2x-if^ , (a;-|)x^ + 2(aix-a2)+ai ^ 
|r(a:: + z7;)| ^ 27re \y\ llH I. 

Using this in (|6.28p . we find that as i? — > oo. 



(6.31) 



d-3 



R^p[RM- Rm, Rm) ^ ^ Li^^-A .. ' 1 ) + ) ' (6-32) 

2"-7T 2 1 j M \ ^ 



where 

^ _ (a;j - i)a;| + 2(aia;-,- - 02) + (.X4 - i)x| + 2(01X4 - 02) + a? 

= S((M + 2m)2 + (M-2m)0"^ ford = 4. (6.33) 

Note that the argument of the exponential in (|6.32p is if M — 2m > 0, otherwise —TrR{2m — M) and, 
in this case, R^p{RM; Rm, Rm) tends rapidly to 0. In all cases, (|6.32p shows that R^p{RM\ Rm, Rm) 

tends to p'^''^kowski(j^^2.^^ ^) (ggg dmni)). 



6.2 Complementary particles 

One benefit of having the explicit formula (|6.25p is being able to examine the case of "complementary" 
particles. The integrand of (|6.27p is meromorphic in n and v. We can rewrite (|6.27p as 

, 5, [ K sinh(7rK) fd—l d—1 d 1 — C 
w.Az, z ) — = ; ; F h in, in : — ; 

^ ^ 7r 2rf+%'^+ir (f ) r (^^) i?2d-4 y 2 '2 '2' 2 

„ f iHi\^ -r-r lEK ieh'\ , d—1 

The integrand is meromorphic in k and i^. It has no singularity when both are real. The Ihs is 
holomorphic in {z^ : v ^ ±i{{d — l)/2 + Z+)}. We analytically continue the integral in the variable 
choose I' complex with Rei^ > and a = Imiy > 0. Recall that, for integer n > 0, 

(-1)" 

z + ^r(z) — (6.35) 

n'.(z + n) 

The poles of the functions k i-^ r(/i ± iK/2) are at k = ±2i{ii + n) {n > integer), and are independent 
of ly. The other poles of the integrand are as follows (n > integer): 

iK fin \ f— 1)" 

+ ^,±iy + n^O ^t(-+P±iv]^— \ > (6.36) 

2 \2 J ^nl{K - 2t(n±iiy + n)) 
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IK f IK \ f— 1~)" 

— +/i±ii/ + n-0 =^r +/i±iz^ — 5—/ ■ -. (6.37) 



2 ' \ 2 ' J -^nliK + 2i{p±iiy + n)) 

The poles k — 2i{fi + iiy + n) = (see ()6.36p ) and the poles K + 2i{ii — iiy+n) = (see ()6.37p ) are on the 
line —2Ileiy + iR. Their mutual distances do not change as v varies, and they all move down as Imz/ 
increases. The poles k + 2i{^ + iv + n) — (see (|6.37p ) and k — 2i{fi — iv + n) = (see (16.36^ ) are the 
opposites of those described before. They lie on 2 Re + iR and move up as Im ly increases. 

If Im v increases from but < Im v < fi, no pole reaches the real axis and the formula (j6.27p continues 
to hold. This is true in particular if fi — ia with < a < (c? — l)/4 = mc/2, corresponding to 
nic > mi > mcV^/2. If this condition is satisfied and itlq > rric, eqs. p.l3p and (|3.14p hold and the 
adiabatic limit exists just as in the case mi > rric- 

When Im v reaches /i we have 

w,y{z, z'Y = \ K p{k, v, v)wk,{z^ z') dK , (6.38) 



Jc 

where the contour C is obtained from R by a small downward excursion to avoid the pole at —2 Re;/, 
and another small upward excursion to avoid the pole at 2Rei^. Once /x < Imi' < /i + 1, we can 
extract the residues of the poles at k = ±2i(/x + iiy). A similar situation occurs when the successive 
poles K = zL2i{fi + iv + n) cross the real axis, so that, for Rei/ > 0, Imv > 0, Imi^ — fi ^ Z, N ~ 
max{j € Z : j < lmi> — /i}. 



Wi,{z, z'Y = I kp{k, z/)w7„(z, z')dK + 

N 

/I . _ I 77 I A' 11/) 

(6.39) 



E 



7, '^2i{p+iv+n)\Zi z)^ W _2i{p+iv+n)\Z : Z) 



Note that if Imi/ < /i [N < 0), the discrete sum is not present. It turns out that AJj(j^) = A„(z^), 
which is consistent with np being even in n. Recall that Wr = w-r for any t. Thus, for N = 
max{j e Z : j < Imv — p}, (always supposing Rei^ > 0), 



N 

W^{z, z'Y ^ / Kp(k, V, l/)Wt,{z, z')dK +^An{iy)w2i(fj.+iu+n){z, z') . (6.40) 

n=0 



We find, for integer n > 0, 



An{v) 



T{2p + 2iv + n)T{~2w - n)T{2p + n)T{-w - n)T{2p + iv + n) 
T{-2n - 2iu - 2n)r{2p + 2iv + 2n)r(i - iv - 7i)r(i + 2p + iv + n) 



(6.41) 



If now we let v tend to ia {Kev tends to 0), (|6.40p will continue to hold provided both parts of the rhs 
remain meaningful. Therefore, if < a < (rf — l)/2, a — /i ^ Z, and N = max{j e Z : j < a — p}, 
p={d~l)/4, 



Wiaiz, z')^ = / Kp{K;ia, ia)w^{z, z') dn + An{ia) 'W2i(a-p,~n){z, z) . (6.42) 

-'^ n=0 



kp{k; ia, ia) = 



Tjp+lf- a)Tip - f - a)np + f + a)r{p - f + a)T{p + ^)r(^ - f ) 

r{tK)T{^iK)T{p + ^ + f)T{p + ^~f) ^ • ' 
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This is obviously positive. For An{ia) we find 

1 r{2a-n)T{2fi + n)r{a-n)r(2n- a + n) 



!2d-%^i?'i-2r(2/i) r(2a - 2/x - 2n)r(i + a - n)r(i + 2fi-a + 7j) 



^ ^ r(2^-2a + 2n) ^ ^ 



The first two factors in this expression are positive since the arguments of aU T functions are positive 
due toa — /X — n>0 and 2fj. — a > 0. The last factor is of the form 

tigi^^(-irn('«)-'^ ("5, 

1 (2n + x) 

^ ' g 

The last product contains n negative factors and the result is positive, so that A„(ia) > 0. Thus the 
Hilbert space with scalar product given by the Ihs of (|6.42p appears as a direct integral of Hilbert spaces 
associated with unitary irreducible representations of Gq. 

We conclude that 

1. Any particle from the principal series can decay into two particles (of equal masses) of any series. 

2. A particle of the complementary series with parameter n — i(3, with < /? < 2/i can decay into two 

particles with parameter ia, a = ^f3 + ^ + where n is any integer such that < n and a < 2fi, 
i.e. n < n — f3/2. This relation can also be written as 

{2^i- (3) ^2(2fi-a) + 2n<2fi. (6.46) 

This implies a form of particle stability, but the new phenomenon is that a particle of this kind 
cannot disintegrate unless the masses of the decay products have certain quantized values. Stability 
for the same range of masses has also been recently found in a completely different context. 



7 Concluding remarks 

In trying to interpret the results concerning the lack of mass subadditivity in the de Sitter universe, 
one can wonder whether they might be due to the thermodynamical properties ([TH [131 IHl [H]) of the 
fundamental state we have been using. We have tested this possibility against a similar computation in 
flat thermal field theory that however does not exhibit this phenomenon in two-particle decays. Another 
issue has to do with energy conservation and the relation mass/energy. dS invariant field theories admit 
ten conserved quantities (in ci = 4). The identification of a conserved energy among these quantities 
has proven to be useful in classical field theory [55]. The same quantity remains exactly conserved 
also at the quantum level although it becomes an operator whose spectrum is not positive [T31 HH [T3] 
even when restricted to the region where the corresponding classical expression is positive [29] : the 
thermodynamical properties of dS fields arise precisely in this restriction [12l [131 [HI [H] ■ Energy is 
conserved also in the decay processes that violate mass subadditivity, once the adiabatic limit has been 
performed. The breakdown of the subadditivity property of masses in dS spacetime just reflects the 
nonexistence of an Abelian translation group and thereby of a linear energy-momentum space. 

When we consider the adiabatic limit problem and its meaning in the de Sitter context a first compli- 
cation is the existence of several choices of cosmic time, having different physical implications and the 
result might depend on one's preferred choice. We have studied the closed and the flat cosmological and 
found that in both models the first factor in (|3.13p diverges like T; thus it has to be divided by T to 
extract a finite result which is the same in both models. 

Here the second (unforeseen) result comes in: in contrast to the Minkowskian case the limiting proba- 
bility per unit of time does not depend on the wavepacket! This result seems to contradict what we see 
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everyday in laboratory experiments, a well known effect of special relativity (Eq. 14. 9p . Furthermore, in 
contrast with the violation of particle stability that is exponentially small in the de Sitter radius, this 
phenomenon does not depend on how small is the cosmological constant. How can we solve this paradox 
and reconcile the result with everyday experience? The point is that the idea of probability per unit 
time (Fermi's golden rule) has no scale-invariant meaning in de Sitter: if we use the limiting probability 
to evaluate amplitudes of processes that take place in a short time we get a grossly wrong result. This 
is in strong disagreement with what happens in the Minkowski case where the limiting probability is 
attained almost immediately (i.e. already for finite T). Therefore to describe what we are really doing 
in a laboratory we should not take the limit T — > oo and rather use the probability per unit of time 
relative to a laboratory consistent scale of time. In that case we will recover all the standard wisdom 
even in presence of a cosmological constant. But, if an unstable particle lives a very long time (>> R) 
and we can accumulate observations then a nonvanishing cosmological constant would radically modify 
the Minkowski result and de Sitter invariant result will emerge. This result should not be shocking: after 
all erasing any inhomogeneity is precisely what the quasi de Sitter phase is supposed to do at the epoch 
of inflation; in the same way, from the viewpoint of an accelerating universe all the long-lived particles 
look as if they were at rest and so their lifetime would not depend on their peculiar motion. 

We thank T. Damour, H. de Vega, M. Gaudin, G. Gibbons, D. Marolf, M. Milgram and V. Pasquier for 
enlightening discussions. U. M. thanks the SPhT and the IHES for hospitality and support. 
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A Appendix. More details in the Minkowski case 

In this appendix we study in more detail the adiabatic hmit in the Minkowski case: it is possible to let 
the two occurences of g in (I3.12p tend to 1 together, or even at different rates. Let 



W(/0, </5l, V2, P)^ j fo{x) fo{y) V2{v) p{<T^) 

Wmo{x, u)wa{u, v)wmo{v, y) dx du dv dy d(7^ . (A.l) 

We will assume that p is C°° and has support in (? + R+, with < c < mo, and that, for each integer 
n > 0, there are constants C„ > and _L„ > 0, such that, for all real t > (? , 

|p(")WI<C„(l + |t|)^". (A.2) 

We take 

^,{x)= f e-'P^^,{p)dp, j = l, 2, (A.3) 

Vj{p) = e-'gjipVej) i^{p°)S{p), ^,{x) = ^ / g,{e,{x^ - t))^ dt, (j = 1, 2), (A.4) 

gj{t)^ [ e-'*'^gj{w)dw, = / e-'*"'V^(w) dw. (A.5) 
Jr Jr 

Here ej = > 0. The function ip belongs to 5(R) with ■(/^(O) = 1. The function g^ is L°° with 

compact support. (The cases of real interest are gj{t) — 0(1/2 — \t\) or gj{t) — d{t)9{l — t).) Wc find, 
after using the various delta-functions, 

^(/O, fl, f2, p) = 

Sip' - m^)0(/) p(^f -pl + mDeih"^) dp dh° . (A.6) 
We now change from the variable h'^ to the variable w such that /i° = + u; : 

e2U{fo, ipi, ip2, p) = e^^ j H{w)gi (^~^^ h (yj^ '^^ ^ J ^^^^^^9i{-r)g2 (^T^^ (^■'') 
Here 

H{w) = (2^)''+3|^;H|2 / {2p'')-^\fo{p)\H{p'-ml)e{p'')p{ml + w{2p'' + w))0{p''+w)dp , (A.8) 



and we set 

Hit) = / e-'''"Hiw)dw . (A.9) 
Jr 

Then 

e2Uifo, Lpi, LP2, p) = {2%)-' Hix)giieix+—y)g2iy)dxdy. (A.IO) 

With our assumptions on p, iJ e 5(R). Since gj is L°° with compact support and H E iS(R), the above 
integral (jA.10|) is absolutely convergent, uniformly in ei and £2- Hence 

£2Z^(/o, ifi, ip2, p) - (27r)-2 /" Hix)Gix, ei, 62) dx, (A.ll) 

G(a;, ei, £3) ^ J 9i (^^i^ + 92iy) '^V ^ j 9i {~y^ 92iy ~ £2x) dy . (A.12) 
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We assume from now on < ei < £2 < 1. Since gj E L°° D and translation is continuous on L^, G is 
continuous in x. 

Example 1. We suppose that gi, 52 are C°° with compact support. In this case the limits when £j 
tend to can be taken under the integral sign in (jA.12p . 

(1.1) if £1 tends to at fixed £2, G tends to the constant gi{0) J 32(2/) dy, independent of £2. 

(1.2) if both £1 and £2 tend to and £i/£2 0, G also tends to gi{0) J 32(2/) dy. 

(1.3) if both £1 and £2 tend to and £i/£2 ^ A e (0, 1], then G tends to the constant / gi{^y)g2{y) dy, 
and 

£2Zi(/o, ^1, p) ^ (27^)""/ H{x)dx J gi{Xy)g2iy)dy. (A.13) 

This holds in particular if £1 and £2 are kept equal so that A = 1. The constant J gi{Xy)g2{y) dy 
may be equal to the preceding constant gi(0) J g2(y)dy, for example if gi{\y) = <?i(0) on the 
support of 92- 

Example 2. We consider the case when gj{x) = 9{x)6{l — x), i.e. gj is the indicator function of [0, 1]. 
Then (see Fig. [2]) 

G(x, £1, £2) = [l + £2x)0{l + £2x)0{-x) + e{x)e{£^^ - e^^ - x) + 

+ £2{£i^ ~ X)e{x - £^^ + £^^)e{£^^ - x) . (A.14) 

G{x, £1, £2) tends to 1 when both £j (with £1 < £2). G also tends to 1 if £1 at fixed £2 and 





1 


















e2' 







-'e-2' 





Figure 2: Graph of G(x, £1, £2) when 51(2;) .92(2;) = 0{x)Q(\ — x). 

then £2 0. Since U e 5(R), the integral (|XTT|) tends to (27r)-2 ^(^.^ g^^, ^ (2txY^H{Q). 

Example 3. Consider now the case gj(x) = 0(1/2— |a;|), i.e. g^ is the indicator function of [—1/2, 1/2]. 
In that case G(x, £1, £2) is even in x (see Fig. [3]): 



G{x, £1, £2) - 0(£r' - £2"' - 21x1) 

£2|x| \ B{2\x\ 



ifl + ^ 
2 V ^1 



£2"^)6'(£j;i +£2"^ - 2|a;|). (A.15) 



G(a;, £1, £2) tends to 1 either if £\ tends to at fixed £2, or if both £j — > (with £\ < £2), and the 
integral (|XTT|) tends to {2?:)-'^ H (x) dx = {2t:)-^H{0). 

Conclusion With the two last choices of gj just described, 



£2U{fo, fi, f2, p) ^ {2^)"+^ p{ml) f {2p°)-'\fo{p)fS{p^ - ml)e{p°)dp . 



(A.16) 
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Figure 3: Graph of G(x, £i, £2) when gi{x) — g2{x) — ^{1/2 — \x\). 

For other choices of gj, the hmit as ei — s- 0, then £2^0 need not be the same as when £1 — £2 ^ 
0. 

If we relax the conditions set on p, the same conclusions hold if e.g. p{s) = p(s; mi, TO2) and d > 4. 

B Appendix. Adiabatic limit (dS): horizontal slices 

Horizontal slices have been described in subsect. 15.11 In this appendix, we study lim^^^+oo T^^Li{fQ, g) 
where g{x) is given by (|5.4p in the coordinates (|5.ip . We denote k = [ttt-q (d- 1)74]^/^ and recall that 
mo > {d— l)/2 (hence n > 0), and Ci{mo, d) = Co(k) (see (I2.27p . Inserting the representation (|2.23p 
for the three occurrences of w,no+ (denoted also in the formula (|3.14p for -Li(/o, g) gives 



1' Co{k) crf,. x,^^ feo(e) K^ite, 9) hoic) dfi^io dfi^ie) 

K.{^,^',g)= f V+.-i (x,0i^~.-..Jx,O9{^)dx. 



ii(/o, 5) 

K{0 = 



(B.l) 

(B.2) 
(B.3) 



We take 7 = 5*0 = G C+ : = 1} ~ S''' ^, the unit sphere in R''. In this appendix, we also set 
i? = 1: a general R can be reinstated in the results by homogeneity. Note that, for any Lp e C^{So x 5*0), 
/ ^{^, V'^d^_j^(a;,OV'Id-i^.^(a;,$')c^M5o(Oc?MSo(C') is in x. Hence, for any bounded g with 
bounded support, if^ is a distribution on 6*0 x So in the variables (^, ^'). We will take g invariant 
under the rotation group in d dimensions (leaving eo invariant), hence Ki^{£^, g) = Kk{L^, L^' , g) for 
every such rotation L. Hence if^ is C°° in ^ when smeared with a test-function in Studying the 
limit of T~^Li(/o, g) for g as in (|5.4p . is therefore equivalent to studying the limit of T~^Ki^{^,^' , g) 
as a distribution in ^' for fixed ^. In this appendix the cases d = 2 and d = A will be treated. The case 
ci = 3, more straightforward than d = 4 (no need to use d = 2), will be omitted. The result in these 
three cases is the same (see (|B.33p and (|B.93P ). 



B.l Case(i = 2 



We use the following parametrizations 



a;*^ ~ sinh t 

= coshi sin6l <( = { ^''^ = - sin (j> (B.4) 

x 



^ — cosh t cos 6 





Particle decay in the de Sitter Universe 



24 



with 

ten, -TT <9 <Tr, {xe X2), -TT < < TT, ^ ^ dV+) . (B.5) 

In these variables, the measure dx takes the form cosht dtd6. For small e > 0, changing t into t ± is 
pushes X into T±. If g{x) ~ go{t), 



JteR. -Tr<e<TT 



iteR, -T:<e<7T 

[smh{t + iO) + cosh(t + iO) cos(6l)]"^"*''[sinh(< - iO) + cosh(i - iO) cos{0 - cosht dtdO. (B.6) 

For real s with < |s| < 7r/2 and real a, 

sinh(i + is) + cosh(t + is) cos(a) = cos(s)(sinht + cosht cos(q!)) + i sin(s) (cosht + sinht cos(q:)) (B.7) 

has a non-zero imaginary part of the same sign as s, so its power n can be taken for any complex fi 
and remains analytic in t for all real t, smooth and periodic with period 27r in a. The integral over in 
(|B.6[) will be performed, using Plancherel's formula, by first computing the discrete Fourier transform, 
in the variable 9, of the two last factors in the integrand, i.e. 

K.i^,e.9) = ^y2 f e'^^^F^t + iO) Fm{t + iO) go{t) cosh{t)dt, (B.8) 

with 

F„,{t + is)= f [smh{t + is) +cosh{t + is) cosi9)]-^-'^ e'^^Ue . (0 < |s| < 7r/2). (B.9) 



By changing 9 to ~9 in the integration, we get: 

F„,{t + ts) = F^^{t + is). (B.IO) 

We use the formula ([H], (15) p.l57) 

P;"(^) = ^l^nil'+l)'^ f ~ ^^"^ cos0]^e™^d0 , {z e Ai, Rez > 0). (B.ll) 

This is stated for m integer and > 0, but by using (see [57], (7) p. 140) 

meZ, (B.12) 



r(A-m+l) r(A + m+l)' 
it is seen to hold for all m € Z. This gives, for < s < 7r/2, 

Fra{t + is) = \2 ^p™^_. (-isinh(t + is)) for t > 0, (B.13) 

r(r7i + ^ — in) 2 

F„,{t + is) = (-l)™e-"/4+^''/2 4^^I^ii— ^P™i_^^(-zsinh(t + zs)) for t < 0, (B.14) 

1 (m + i — ik) 2 

Fm{t + is) = e-/4+-''/2-?^^^^±^p-™ (*sinh(t-zs)) for t > 0, (B.15) 

1 (— m + i + zk) -2+^1^ 

= (-l)"e"/^+-/^ ^^^I , P''\^J'^^'Mt ^s)) fort<0. (B.16) 
Taking go{t) = 6l(T/2 - t)9{t + T/2), we rewrite (iRSl) as 

E^""'(^™+^™)' (B.17) 
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where 



/+ = 7^ / ^ Emit + iO) Frn{t + lO) COSh{t)dt 

Jo 



27re''''(-l)'" 



sinh(T/2) 



(B.18) 



(B.19) 



In the last expression we have used T{^ + z)T{^ — z) — it / cos(7rz), and changed to the variable u = sinhi. 
Similarly (with now u = — sinh t) , 



1 /-o 

/- ^ — / Fra{t + iQ) F^{t + COs\i{t) dt 

J -T/2 



<.sinh(T/2j 



(B.20) 
(B.21) 



Remark B.l Using (ah. 10) and r(i + z)r(i ~ z) = tt/ cos{ttz) shows that the rhs of (|B.21[) can be 
obtained from the rhs of (jB.19|l by changing in the integrand (but not outside the integral) k to —k. 
Note also that /± = by (|B.10p . 

The meaning of (|B.17|) is that KK,{^,^',g) is a distribution in ^' as expressed in the coordinate (j), and 
that m 1-^ + is its discrete Fourier transform. It is tempered, i.e. |J+ + /~| does not increase 
faster than a power of \m\ as \m\ oo. To prove that T^^K^ tends to a limit (also a distribution in </>) 
as T — > oo is equivalent to proving that 

(1) For each m, T^^(/+ + Z,7j) tends to a hmit Um as T ^ oo, 

(2) there are two positive constants P and Q such that T^^|/+ + /~| < P(l + |m|'3) for all m and T. 
If both conditions are satisfied, Um is the m'^ Fourier coefficient of the limit, i.e. limT^^i^K = 



B.1.1 Condition (1) 

We need the asymptotic behavior of P'^{z) as \z\ oo, as described in ^57], pp. 123, 124, 126, and 164. 
For z G Ai and C = z~^. 



2-A-l^-l/2p(_ 1 _ ^) ^-A-l+M (^2 _ ^)-p/2 

r(-A-M) 

2^7r-V^r(i + A) (z^ - 1)-^/^ 

r(i + A-^) 



P(i + A/2 - fi/2, 1 + A/2 - /i/2; A + 3/2; C) 



P(-A/2 - ^/2, 1 - A/2 - m/2; 1 - A; 0- (B.22) 



If A + i ^ Z, the two hypergeometric functions can be expanded into convergent power series for |^| < 1. 
For z e Ai and \z\ oo, we find 



2-A-i^-i/2p(_ 1 _ ^) 2V-i/2r(i + A) Z^ 



r(-A-M) ' r(i + A-M) 

Hence, as u — s- +oo, 

2-i+iK7r-i/2r(jK) e*^/4+'^«/2 2"^"'''7r"i/2r(-iK) e*'^/4-'^«/2 u-^-' 



(B.23) 



r(i + iK — m) 



r(i — in — to) 



2-|-iK^-l/2pj^_^^-j g-i7r/4+7rK/2 ^-i-iK 2~2~'"'''7r~"'"/^r(iK) g-*'^/4-irK/2 y-j+it 



r(i — iK + m) 



r(i + in + to) 



,(B.24) 
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We first consider the 'ofF-diagonal terms' of this product: 



r(i + in — TO)r(i + in + m) — in — m)T{^ — in + m) 



(B.25) 



These two terms are exchanged by changing k to — k. The contribution of the first to I^-^/T is of the 
form 

, .sinh(T/2) , 

Const.-/ u-i+2"'dM = Const.— — (sinh(T/2)2'«-l). (B.26) 

This tends to zero as T — > +00. The same happens for the second term. The 'diagonal terms' are 



r(i + IK — m)r(i — IK + m) r(i — in — m)r(i + + m) 

Again these two terms are exchanged by changing k to — k. Their sum can be reexpressed as 

(--l)'"cosh(7rK)2M-i 
7r«;sinh(7rK) 

Since 

/.sinh(T/2) 



(B.27) 



(B.28) 



/snin^_i / z ) 
u-^du = log(sinh(T/2)) - T/2, (B.29) 

J 7rKsmh(7rKj 



Because of Remark lB.ll I„-^/T has the same limit as I^/T and 

1 , , , 26^*" cosh(7rK)2 

[/„ = lim - /+ + = . ■ / , . B.31 

T-.+00 r K smh(7rK) 

C/m is independent of to, so that if Condition (2) is satisfied, 

1 ,^ ^1 s 47re''"cosh(7rK)2 ^, 47r e'"' coshfTrK)^ ^ _ 

hm - K^it r , 5t) = . , / / 5(0) = . , / / it e) , (B.32) 

T— ►+00 i K smh(7rKj K sinh(7rKj 



(B.33) 



and (see (|B.ip 

r rr-ir tf \ 2r'f \ 47re^'' cosh(7rK)^ 7^ tt coth(7rK)^ 

Inn T ii(/o, gr) = 7 C'o(k)c2,k r-rT — ^ = n ' 

K smh(7rK) \k\ 



We note that, owing to the delta function in (IB.32|) the dependence on ho (i.e. on /o) has completely 
disappeared from the limit. This result agrees with (|5.9|) . 

B.1.2 Condition (2) 

In this subsubsection A always denotes — in with k S R and k 0. We first return to the first 
step of the preceding subsection in the case to = 0. From the identity (jB.22p and the analyticity of 
1-^ F{a, b, c; Q in the unit disk, it follows that there is a Mo{k) > such that 



< Mn{K)\z\-'^/^ zeAi, |z|>2. 

(B.34) 



By (IB.13IIB.14|) . there is also an Mi(k) > such that, for < s < 7r/2, 

iFoit + is)\ < Mi(k) I sinh(f + is)\-'^^^ . (B.35) 



Particle decay in the de Sitter Universe 



27 



We now obtain crude bounds for |F„j|. For t > and < \s\ < 7r/2, changing 9 to 6 + tt in (|B.9[) . we 

get 



F,n{t + is) = (-1)™ cosh(i + is)^ / ((1 - cos(6l)) - (1 - th{t + is))f e™^ d9. 

J — 7r 

Changing to the variable ip — 9/2, 

Fm{t + is) = (-l)™2(2cosh(i + is))^A™(z), 

_ , (l-tht)(l-2tgs) 

1 + i th i tg s 



= i(l-th(t + is)) 



We now suppose < tgs < 1/4. It foUows, after some calculations: 



tgs < 



Im z 



Rez2 

We define z = x — iy with x > Q. Then 



< tg(2s), tg(s/2) < 



Im z 



Re z 



<tgs, 



(B.36) 

(B.37) 
(B.38) 

(B.39) 
(B.40) 



< a; < |z| < 3/4, < y < a; tg s < 3/16, — ^ < \z\ < x J I + ig^ s < x^/YJ/TQ. (B.41) 

V2 



Recall that for p > 0, -tt < 61 < tt, C G C, 
Thus 

e-"l"l|^„(z)| <il(z) = / |sinV-2'r'/'d<y5 = 2 / 

Jo 



dt 



(B.42) 



(B.43) 



After splitting the integration interval as [0, 1] = [0, a;]U[a:, V3/2]U[\/3/2, 1] , straightforward estimates 
give 

Aw _ 

(B.44) 



e~''l''U„(z) < H{z) < ^= + 8 + 41og(V3/x) < 27 + 4i 
3v 3 



We now consider 



Ao(z)-A™(z)= / (sin2((p)-z2)A(^_g2»mv)^^ ^ / (sin2(^)-z2)^(l-e2'"^+2imv3)d(^. (B.45) 



Using 



we get 



|1 - e^™'^ + 2i'mLp\ < 2mV^ < -^sm^{ip) Wip e [-tt/2, tt/2] 



(B.46) 



-ttIkIi 



2 2 




TO TT 








2 2 




TO^TT^ 




2 










z2 



r/2 



.2|-l/2 «^„2/ 



7r/2 



< 



iJ(z) + TOVM |sin2((^) - 22|l/2 



(B.47) 
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Since |z^| < 17a;^/16 and 2x^ log(l/x) < 1/e, there is a constant AI2 > such that, for aU rn, 



{B.4 



and hence 



\Fo{t + ts) - (-l)"i^„(t + is)| < V2e2''l''l|cosh(t + is)-i/2-«|^,^^^2 ^g ^g^ 

With < tgs < 1/4, as we have chosen, | cosh(t + is)-i/2-*«| < (cosht)-i/2(17/16)i/'^el''l/'^, so 



that 



\Fo{t + is) - {-irF,n{t + is)\ < V2e(2-+i/4)l-l(i7/i6)i/4|cogl^(^)|-i/2^^^^2^ (g_5Q) 
and, by (|B.35p . there is an M3{k) > such that 

\F.m{t + is)\ < M3{k){1 + m^)! sinh(t)ri/^ (B.51) 



Therefore, using the bound (|B.43[) . independent of to, for < t < ti, and the bound (|B.5ip for 
ti < t < T/2, we find that there is a constant M4(k) > such that 



IT-I/+I < Af|(K)(l +to2)2 VtogZ. 
The same holds for \T~^I~-^\. This proves that Condition (2) is satisfied. 



(B.52) 



B.2 Other dimensions 

In this subsection the dimension of the de Sitter space-time X is n — d > 2, i.e. the ambient Minkowski 
space-time is R'^+i. The notation n = d is used to stay close to Chap. IX, p 448 ff, which is 
constantly used in this section. As in Subsect. IB. 1.11 we wish to compute 



(B.53) 



X 



where ^, C' G dV+ C R'^+^ and ^" = = 1, A = -{n - l)/2 - iz/, ^f{x, are as defined in (1^151) . We 
use the following parametrization for x — x) £ X, ^ = (1, ^) e ^' = (1, G 9V+ (see [3D], 

p. 448]). 





x*' = sinht, a; = — coshtu 






= 1 




' = 1 






= sin0„_i . . . sin 02 sin0 


1 


e 


= 




= sin(?!)„_i . 


. sin (j)2 sin t/fi 




u'^ ~ sin 6n-i ■ ■ ■ sin 62 cos 6^ 


1 








^'^ = sin (/)„_! . 


. sin (l>2 cos 01 


< 


M'^ — sin0„_i . . . cos 6*2 


< 






< 


^'^ = sin(/)„_i . 


. cos 02 




= sin0„_i cos0„_2 
w" = cos6'„_i 




e 
[ e 


-1 = 
= 1 




^ =sm0„_ 
4 = cos(/)„_i 


1 COS0„_2 


Here i € R, < 6*1 < 27r, < 


< 27r, < 


Ok < n 


for fc > 1, < 0fc < 


TT for fc > 1. 



(B.54) 



notations 

dx — cosh"^^ t dt du, du = sin"^^ 0„_i d0n-i ... sin 02 d92 dOi 
We also use the normalized measure da{u) on S"~^, 

27r"/2 

dii — 

/S" 



da{u) — du, r2„ — du — 



r(n/2) • 



(B.55) 



(B.56) 



We restrict g to be of the form g{x) = grix) = go{t) = e{T/2 - t)e{t + T/2), T > 0. The integral ((B331) 
takes the form 

K^{^,S,',g)^ I go{t) {coshty-^ dt [ F{t, u)G{t, u) du , (B.57) 



R 
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F{t, u) = (sinh(< + iQ) + cosh(i + iQ) cos(6'„_i))^ . (B.58) 
For G, we have G{t, u) = {x+ ■ ^')^. Note that ^' — R^, where R is the rotation in R" 

For example 



/ 1 



M„„_i 











\ 



(B.60) 



... 1 

... cos0„_i sin(/)„_i 
\0 ... -sin(/i„_i cos(?;>„_i / 

Therefore 

G(t, u) = F{t, R~^u). 

As in the case d = 2, we reexpress the integral over u in (|B.57|) using harmonic analysis on the 
sphere. 

Harmonic analysis on 5"^^ uses an orthonormal basis {S^^-} of functions on the sphere = 0, 1, 2, . . ., 
if is a multiindex) . This is fully described in [30J, Chap IX]): 



(B.61) 



(B.62) 



For fixed i the functions {S^} generate a finite-dimensional subspace 7i"^ of L^(5'"~"'^) in which the 
regular representation of SO(n) reduces to an irreducible unitary representation, characterized by its 
matrix elements in the basis {S^}: for any g S SO{n), 



(B.63) 



M 



The functions {^^} and {t^KAi} ^^'^ analytic on S"' ^ and SO{n) respectively. Given two arbitrary 
functions hi, /i2 on S"""^ we have (for j — \, 2) 



(., K 



hj (u) S]^ (u) dij(u) , 



/ii (m) /i2 (u) da{u) — 



/ill- 



These formulae imply that 



where (55,.-i(u, u) denotes the distribution (actually measure) on S"" x S" defined by 



5511-1 (u, v)ip{u, v)dudv 



S"-ixS"-i 



(p{u, u) du 



(B.64) 

(B.65) 
(B.66) 

(B.67) 
(B.68) 



Actually, as is the case for all invariant distributions on 5" ^ x 5" ^, smearing 55,i-i(-u, u) only in w 
with a C°° function produces a C°° function of u: 



Sgn-i (u, v) ijj{v) dv — il^{u) 



(B.69) 
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Choosing in particular u — e^, we can use the formula ([30], IX 4.1 (1-4) and text there) 



Inserting this in (|B.67p gives 



lT{e + n~2){2e + n~2) „^ ^ 

eiT{n-: 

Taking v = with f' given by (|B341) and using [30j IX 4.1 (3-4) 



/ -^ T{e + n-2){2£ + n-2) „ 

Slf (5e„) = iiT(n_l) ^^'^^^ 

we get 

A ^-'^ i^r(^+n-2)(2£ + n-2) , 

with R given by (jB.59|) . Harmonic analysis extends to distributions on the sphere, as it does on S^. 
We apply (jB.64IIB.66|) to the case hi{u) — F{t, u), h2{u) = G{t, u). Because F{t, u) depends only on 

cos On-l, 

fK{t)= I F{t, u)~'^{u)do{u)^5Kof,{t) . (B.74) 

Note that t can be complexified in (IB.74p . i.e. t can be replaced by t+is with < |s| < 7r/2. In the sequel 
we omit the i-dependence of /^(t), writing simply unless the t-dependence becomes significant. We 
have 

G{t, u)^Y.fo^oiR-'u) = J2{f^tUR))^Kiu) ■ (B.76) 

£ £, K 



Therefore 



Also 



F{t, u)G{t, u)du = ^4(i?)|/oT, (B.77) 

-1 

K.ii,e,9T) = n^J2^UR) \foit)\'icoshtr'Ut (B.78) 

e J-T/2 

For (IR751) see [30] IX 3.6 (6,7) p. 480. For (lB79l) see (S^, IX 4.2 (8) p. 484. We thus have 

/o = ^n^Ao j [sinh(i + iO) + cosh(t + iO) cos6i„_i]^ C7^(cos6'„_i) sin""^ 6»„_i dOn-i ■ ■ ■ sin 02 d92 d9i 

= n^^n^-iA^o / [sinh(t + iO) + cosh(t + iO) cos 0]^ sin""^ ^ (cos 0) d0 (B.80) 
Jo 

In these formulae Cg is a Gegenbauer polynomial: see [57] p. 175 for the definition. The formulae [57] 
p. 176 (9), and [30], IX 3.1 (3), giving the explicit coefficients of coincide, so we are dealing with 
the same objects. 



Particle decay in the de Sitter Universe 31 

B.3 The case d = n = 4 

Wc now restrict our attention to the case d = A, keeping the notations of the preceding subsection. In 
this case A = — 3/2 — jk, O4 = 27r^, Aq = 1. We exclude the case k — 0. Since {n — 2)/2 = 1, the formula 
(|B.80|) gives: 

f^^- [ [sinh(t + iO) + cosh(t + iO) cos 6]^ C} (cos 9) sir? 9 d9. (B.81) 
7^ Jo 

We have ([57], 3.15.1 (15) p. 177) 

Cl(cos^) = 5^^fc^. (B.82) 
smd 

Therefore, for sufficiently small s > 0, 
2 

f^{t + is) = - [smh{t + is) + cosh{t + is) cos 9]^ sm{i + 1)9 sin 9 d9 (B.83) 

Jo 

1 

= -cosh{t + is)^ [th(t + is) + cos6']^ sin(^+ Ije* sin6'(i6' (B.84) 

I" J-TT 

- ^ '^^ cosh{t + is)^ [ [thit + is) + cos 9]^+^^ cos{£ + 1)9 d9 (B.85) 



7r(A + l) 

= '•f cosh(t + is)-'^ [ [sinh(t + is) + cosh(i + is) cos 9]^+^ cos(£ + 1)9 d9BM) 
7r(A + l) 

Recall that A + 1 = -1/2-iK. Therefore, comparing (|B.86p with (|B.9|) . we find, for < |s| < 7r/2, 

+ - 2.(-^ ^^l|cosh(. + ^s) + + ^-(^+^)(^ + 

^^^^^ —F,+,{t + ts), (B.87) 



7r(— i — zk) cosh(f + is) 
using ((Biol) . Hence 

-r^ (£+1)'^ f'^^^ 

K^tC^gr) = n^Y^tUR) 2 2^ /.A / ^>+i(i + iO)F,+i(i + zO) coshidt (B.88) 

= E4(i?)^^^^Uii+/r+i) (B.89) 

with the notations of (jB.18p . Therefore, by (|B.3ip and the proof of Condition (2) for d = 2 (subsect. 
IBTD . 

J (k^ + 1/4)k smh(7rK) 

In the case n = 4, (jB.73|) becomes 

^(/ + 1)2 4(i?) - n^Ss^iei, e) = in'Ss^iC e) • (B.91) 

£ 

Thus 

1 , , , leTT^e'^" cosh(7rK)2 , , , 

T^+oa 1 (k^ + 1/4)k smh(7rK) 

It follows (see (jB.ip that 

r rr.-ir fr ^ 727rcoth(7rK)2 

hm T 'Liifo, gr) = n • ^-^^^ 

T-t+oo \k\ 

This is the same as in the case d — 2. 
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C Appendix. Adiabatic limit (dS): parabolic slices 



We again take R ^ 1. We again start from the formulae (|B.mB.3[l of Appendix [B] but we only require 
K € R \ {0}. The function g will be chosen as announced in subsect l5.ll The map {t, y) i— s- x{t, y) 
defined in (|5.5p is a diffeomorphism of R'' onto the "upper half" X^'^ = {x G Xd : + x"^ > 0}, 
and (t, y) i— > —x{t, y) is a diffeomorphism of R'' onto the "lower half" = —X^^. The cycle 7 

appearing in (jB.ip will be chosen as Vb = C+ n G Md+i ; C° + = !}■ It can be parametrized by 
the diffeomorphism 77 1-^ ^(77) of R'^^^ onto Vq : 




^1 + ^/'), 

{l<j<d~l). 



(C.l) 



Thus Vo is a Euclidean space with (d^ • c?^) — —dtf' on Vq- The stability group of the vector ep — in 
Go leaves Vq invariant and acts as the group of Euclidean displacements there. As noted in Remarks 12. II 
and l2.21 the Gq invariance and homogeneity of il)^{x, ^) imply that it can be regarded as a distribution 
in ^ on Vb, C°° in x on Xd- For a real g G iS(Xd), if we denote g{x) = g{—x), we find 

K.{^, 9) = e'^-K^^{^,^',g) . (C.2) 

It will turn out that g can be chosen invariant under the stability group of eo — e^. Then K,^ is an 
invariant distribution on Vq xVq. For our purposes it will suffice (and be possible) to study the limit of 
T-1X,(^,C',<?t) with 

grix) ^e{t + T/2)9{T/2~t), t = log(a;° + x'^) , (C.3) 



and to add in the end the limit of T ^Ki^{£^,£^' , gx) obtained from (jC.2[) . 
With X parametrized as in (|5.5p and ^ as in (jC.ip . we have 

(y - ,7)2 o 



For k e R'^ \ we find 



2s 



— s = e 



(C.4) 



7^^„ (k, 77) = / y), Oe^'^dy 



iky , 



d-l - -, 

2— iv^zkr] 



2 



(s T ie) 



dy 



a 



(s T «e) 



= 2° 



where A; ~ |k| and s = 



-1 — 21/, 



2/ 



R<i-i 



dy 



(s T «e) 



y 2 Ji^{ky)dy , (C.5) 



We use the following formula ([31j (51) p. 95 with some notational 



changes) 



(y 



2 ' z^) 



-^+''^V^Jx-i{ky)dy^ [ '- 



(C.6) 



/Q y ^ / r(A — ia) 

This is valid provided fc > 0, Re^; > 0, ReA > 0, and Re(A — 27a + 1/2) > 0. Note that none of these 
parameters except k needs to be real. In our application, A = (rf — l)/2 and a — v £ \ {0}. We 
take z = —is + e, s > 0, e > arbitrarily small, a = v. Since {y"^ + z^) then has a small negative 

imaginary part, this will correspond to the case of "(Ma+w (|C.5|) . In (jC.6|) K^ia denotes the Macdonald 
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function. This will introduce no lasting ambiguity since we will use the identities ([31], (5), (6) p. 4, 
(15) p. 5) 

K_U-iks) = -^e^^^'^H^Xiks) = ^ , (J..(fcs) - e-^'^J.^ks)) . (C.7) 

2 2sinh(7ri/) 

This yields 

(k, = e-^-' ^g^;^!^;^^ (e^^J„(fc.) - J_(fc.)) , A . . (C.8) 

■!/']'^jj^_jj^(k, s, 77) can be obtained from this since, for real t^, it is the complex conjugate of i/jZA+iiyC^k, s, 77): 
,A:,_..(k, ,) ^ e-^'^ ^^^^;^^^^^^^^^^^ (e- J_.(..) - J.Aks)) , A _ . (C.9) 

Remark C.l By the preceding remarks, if ^ € Vq is expressed in terms of 77 e R"*"^ as in (|C.1[) . ^^{x, ^) 
is a tempered distribution in r/, a C°° function of x, and an entire function in a. If x is expressed as in 
(15. 5p , its Fourier transform with respect to the variable y is also a tempered distribution in the variable 
k conjugated to y and in 77, C°° in s and holomorphic in a, and, in this sense, the formulae (jC.8[) and 
(|C.9[) can be continued to all v. If v is taken real in these formulae, their rhs becomes locally bounded 
in k, in particular locally . 

Supposing g{x) = G{{x^ + x'^y^) (for example if g{x) = grix) Gt{s) = 9{s - e~'^ l'^)Q{e^ ^"^ - s)), 
Plancherel's formula gives 

K.{i,i:,9) = (^i.f-" f s-^G{s)i;~ (k, s, (-k, s, 7j')dsdk . (C.IO) 

Inserting (jUTS)) . we find 

K,{^,e,9)^i2^y-' [ e*-(''-''')Fjk, 5)dk, (C.ll) 

/■OO 

i^Jk, 5) = / s-iG(s) J,,(fc5)J_,,(A:5) + BJliks) + CJl,,{ks)] ds , (C.12) 
Jo 

where 

^ ^ 2'^7r'^+ie""cosh(7r^) ^ 

r (^ + iK) r - m) sinh2(7rK) ' 



r (^ + »k) r (^ -l^sinh^(7rK) 



B^C^—r-, ^ — -r-^ ^ . . (C.14) 



Going back to eq. (|C.12|) . we divide the integration range into the intervals [0, 1] and [1, 00]. After 
dividing by T, the contribution of the second interval is bounded by 

Const.-/ fc^^s^^ds = Const. — . (C.15) 
T J I kT 

This is because | jQ(a;)| < Const. a;~^^^ as a: ^ +00 (see [31], p. 85). Hence the contribution of the 
second interval tends to as T tends to +00. The function Jq, can be written as 

Mz) = (z/2)" f; = (,/2)« + ■ (^-i^) 

^-^ m!i mH-a + l) \i(l + Q;) / 

m— ^ ^ \ \ / / 
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Thus as T tends to 



-T/2 



e-T/2 
A 



[AJ,,{ks)J^Uks)+BJl{ks) + CJ'i,^{ks)\ ds 



B{ksl2Y"' C{ks/2) 



-2iK 



Const. A:^s^ 



ds 



A 



T{l+iK)T{l-iK) r(i + iK)2 r(i-iK)2 

Sfc2«(l - e~^^'') C/c-2«(l-e*^'=) Const. P(i _ e-^) 



2r(i + m)r(i-iK) 2irKr(i + iK)2 -2irKr(i- iK)^ 



2r 



Hence 



This gives 
and (see (IRT|) ) 



1 ~ ^ ^sinh(7rK) 

rli^oo T^'^^^' = 2r(l + z^)r(l-z«) = 27r^ 



1 N Asinh(7rK) ^, ,, 
hm ^^^S{v-i/), 



hni —Li ( fn, Or) = 



7^Co (K)c£i_K A sinh(7rK) 7^7r coth^ (^^) 



27rK 



2|k| 



(C.17) 
(C.18) 
(C.19) 

(C.20) 



This is half of the result in (jB.33p or (|B.93p , but it is doubled by the addition of the contribution of 

9T- 



D Proof of the projector identity 

In this appendix, we give a proof of the formula (|2.25p in the de Sitter case, with masses m and to' in 
the principal series, i.e. w? = ^"^ + {d - 1)^/4, to'^ = v"^ + {d - 1)^/4, with real jJL ^ Q and y ^ Q. We 
set R equal to 1. The meaning of (|2.25p is 

lim / W,n{z, x)W^,{x, y)g{x)dx^Ci{m, d)5{m^ -m''^)Wrn{z, y). (D.l) 

For g G S{Xd) the integral in this formula is well defined (see Remark 12. 2p . The same method as in 
Appendix [Cl will be used. Using Eq. (|2.23p reduces the problem to the study, as g tends to 1, of 

K^Ate,9)= f Ax.0r.-^^Ax.O9{^)dx . (D.2) 

Recalling Remarks 12. 1112. 2 1 and IC.li and parametrizing ^ and ^' in terms of 77 and rj' as in (jC.ip . we see 
that, for a general smooth fast decreasing g, this is well defined as a distribution in 77 and 77', and an 
entire function in ^ and u, and, denoting g{x) — g{~-x), it satisfies 

K^.A^, e, g) = e-(A'+-)if_^_^,(^,^',g) . (D.3) 

(It is sufficient to verify this formula for real /i and v). We will use the same coordinates (|5.5p and 
many of the formulae of Appendix [Cl We wish to take g as gu{x) = 0{x^ + x"^), 01 gd = gu- Thus gu {u 
stands for "upper") is the indicator function of the domain covered by the coordinates (|5.5p . We denote 
KJIA^'^') = K^A^^i'-,9u) and Kf^A^^^') = ^t^Ai^^' ^9u)- To make the integral converge, we first 
replace gu by a better behaved g^ of the form g^{x{t, y)) = Ge{e~*)guix) which will tend to guix) as 
£ ^ 0. We thus consider 

KA^.i')^K^Ai.i'.9l)^ j^^^ ^Pt^_^^^^{x,0i'Zi^+J^^C)9Ux)dx , (D.4) 
KtA^^a - K^A^,e,9l) = e^<^^^+''^Kl^l_,i^,e) . (D.5) 
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We now take /i and v real and furthermore require ^lv > 0. Using the coordinates (j5.5p and parametrizing 
^ and ^' as in Appendix [Cl (see (|C.1|) ). we may use the Plancherel formula as was done there. We 
obtain 



(D.6) 



-id-l 



■ X 



sinh(7r^) sinh(7rz^) T -ifj)T + i j/) 

■s 



(D.8) 



where k = |k|. We can use the following formula ([31j, 7.7.4 (30), p. 51): 

Ja{as)Ji3{as) s^'^ ds ~ 

(a/2)''-ir(p)r((a + /3 + l-p)/2) 



2r((l + a + f3 + p)/2)r((l -a + P + p)/2)r((l + a-(3 + p)/2) ' 
Re(a + /3 + l)>Rep>0, a>0 . 



(D.9) 



Choosing Ge(s) = with < e < 1, and using (|D.7p to obtain i4r^;J(k) from i4r^:J(k), we obtain 

^rk^ + P^rk^ - 2''-V+^fc'(^-)(fc/2)-^r(i-£) 

(e'^(M+'^) + l)r((-i;i + ijy + e)/2) 
r((2 + e)/2)r((2 + i/z + w - e)/2)r((2 -i^i-iv- e)/2) 
(gTTM _^ e'^'^)r((-i^ - w + £)/2) 



r((2 — ip, — iu 


-e)/2)r((2- 


'r ijJL — iv - 


-e)/2)r((2- 


- i/i 


+ iiy 


-£)/2) 






ii^ + e)/2) 








r((2 + + ii^ 


-e)/2)r((2- 


- ifi + ii/ - 


-e)/2)r((2 4 


- 


— iv 


-£)/2) 






-i)r((z/i- 


- w + e)/2) 








r((2 + ifj, — iv 


-e)/2)r((2 


— — iv - 


-e)/2)r((2- 




+ w 


-£)/2) 


have well-defined limits in 


the sense 


of distributions 


in /i 


and V. 



(D.IO) 



of each term inside the square brackets we make the substitution r(z) = r(l + z)/z. As e — + 0, we 
find 



sinh(7r/i) sinh(7r:/) T - i ii)T + iv) 



r(i + iii±^)r(i-iii± 



1 1 



r (1 + i^) r (1 - ^^^2^) + i^) + e -i(Ai + i^) + e. 
Using [it + £)""'" + (—it + e)^"'" ~ 27r (5(i), and r(l + iz)Y{\ — iz) = ttz/ sinh(7rz), this gives 



(D.ll) 



Particle decay in the de Sitter Universe 



36 



d[fi — I'j -\ d[fi + f) 



(D.12) 



Recall that we are interested in the case when /i 7^ and ^ have the same sign. In this case 
(5(/i + i^) = 0, and \fi\^^6{fi — i^) = 2S{fi^ — v^). Thus, in this case, 



(D.13) 



Therefore, by jES]) and (|DJ)) . 



2rf+2^<i+ig-M|coth(/x)| ^, 2 /^ 



(D.14) 



Recall that this holds when [i and v are both non-zero and have the same sign. Still in the same case, 
using (I2.23P we have 



x')dv 



(D.15) 



7x7 



We choose 7 = Vo as described at the beginning of Appendix [C] and of this Appendix. With the 
parametrization (|C.1[) . this is a (d— 1)-Euclidean space and dn-yiO = d'^^^i]. Therefore, by (|m4)) . the 
rhs of (|D.15P is given by 



2<i+2^<i+lg7rM|coth(/i)| 



and finally 



W^ix, y)>V.(y, x') dy ^ 27r\coth{fi)\S{ti^ - iy^)W^{x, x') . 



(D.17) 



Xa 



Although we have assumed /i and v to have the same sign in the derivation, it follows from = 
and the form of the formula above that it holds for all possible relative signs, provided fJ. ^ and 
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